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, There is now firm evidence that the ICM consists of a mixture of hot plasma, mag- 

' netic fields and relativistic particles. The most important evidences for non-thermal 

phenomena in galaxy clusters comes from the spectacular synchrotron radio emis- 
sion diffused over Mpc scale observed in a growing number of massive clusters and, 
more recently, in the hard X-ray tails detected in a few cases in excess of the ther- 
\ mal bremmstrahlung spectrum. A promising possibility to explain giant radio halos is 

given by the presence of relativistic electrons reaccelerated by some kind of turbulence 
generated in the cluster volume during merger events. With the aim to investigate the 
connection between thermal and non-thermal properties of the ICM, in this paper we 
develope a statistical magneto-turbulent model which describes in a self-consistent way 
the evolution of the thermal ICM and that of the non-thermal emission from clusters. 
. Making use of the extended Press & Schechter formalism, we follow cluster mergers 

and estimate the injection rate of the fluid turbulence generated during these energetic 
events. We then calculate the evolution of the spectrum of the relativistic electrons in 
the ICM during the cluster life by taking into account both the electron-acceleration 
O ' due to the merger-driven turbulence and the relevant energy losses of the electrons. 

We end up with a synthetic population of galaxy clusters for which the evolution of 
c/3 ■ the ICM and of the non-thermal spectrum emitted by the accelerated electrons is cal- 

^ ' culated. The generation of detectable non-thermal radio and hard X-ray emission in 

the simulated clusters is found to be possible during major merger events for reliable 
, y values of the model parameters. In addition the occurrence of radio halos as a function 

rN I of the mass of the parent clusters is calculated and compared with observations. In 

^ . this case it is found that the model expectations are in good agreement with observa- 

tions: radio halos are found in about 30% of the more massive clusters in our synthetic 
population {M>1.8 x lO'^^ M©) and in about 4% of the intermediate massive clusters 
(9 X 10^"* < M < 1.8 X IO^^Mq), while the radio halo phenomenon is found to be 
extremely rare in the case of the smaller clusters. 

Key words: acceleration of particles - turbulence - radiation mechanisms: non- 
thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

1.1 Evidence for non-thermal phenomena 

There is now firm evidence that the ICM is a mixture of hot 
gas, magnetic fields and relativistic particles. 

The most important evidence for relativistic electrons in 
clusters of galaxies comes from the diffuse synchrotron radio 
emission observed in a growing number of massive clusters. 
The diffuse emissions are referred to as radio halos and/or 
radio mini-halos when they appear confined to the center of 



the cluster, while they are called relics when they are found 
in the cluster periphery (e.g., Feretti, 2003). 

Diffuse radio emission is not the only evidence of non- 
thermal activity in the ICM. Additional evidence, comes 
from the detection of hard X-ray (HXR) excess emission 
in the case of the Coma cluster, A2256 and possibly A754 
(Fusco-Femiano et al., 1999, 2000,03,04; Rephaeh et al., 
1999; Rephaeli & Gruber, 2002,03). If these excesses are in- 
deed of non-thermal origin, they may be explained in terms 
of IC scattering of relativistic electrons off the photons of 
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the cosmic microwave background (CMB) (Fusco-Femiano 
et al., 1999,2000; Rephaeli et al., 1999; Volk & Atoyan 1999; 
Brunetti et al.2001a; Petrosian 2001; Fujita & Sarazin 2001; 
Blasi 2001). Unfortunately the poor sensitivity of the present 
and past HXR-facihtics docs not allow to obtain a iron-clad 
detection of these excesses and thus future observatories 
(e.g. ASTR0-E2, NEXT) are necessary to definitely con- 
firm them (see Rossetti & Molendi 2004; Fusco-Fcmiano et 
al.2004), to constrain their spectral shape, and hopefully to 
increase the statistics of HXR excesses. 

The presence of high energy hadrons is at present not 
yet proven, but in principle, due to confinement of cosmic 
rays over cosmological time scales (e.g., Volk et al. 1996; 
Berezinsky, Blasi & Ptuskin 1997), the hadrori content of 
the intracluster medium might be appreciable and may be 
constrained by future gamma-ray observations (e.g., Blasi 
2003; Miniati 2003). 

1.2 Thermal — non-thermal connection 

Giovannini, Tordi and Feretti (1999) found that ~5% of clus- 
ters in a complete X-ray flux limited sample (from Ebeling et 
al.l996) have a radio halo source. The detection rate of radio 
halos shows an abrupt increase with increasing the X-ray lu- 
minosity of the host clusters. Indeed it has been found that 
about 30-35% of the galaxy clusters with X-ray luminos- 
ity larger than lO**^ erg s~^ show difi'use non-thermal radio 
emission (Giovannini & Feretti 2002). The high luminosity 
of clusters hosting radio halos implies that these clusters 
also have a high temperature {kT >7 keV) and a large mass 
(^lO"''' Mq). Although still based on relatively poor statis- 
tics, these observations suggest a leading role of the cluster 
mass and/or temperature in the formation of radio halos. 

In addition correlation between the radio power of radio 
halos at 1.4 GHz and the bolometric X-ray luminosity, tem- 
perature and mass of the parent clusters have been found 
(Colafrancesco 1999; Liang et al. 2000; Giovannini & Feretti 
2002; Govoni et al.2001). 

1.3 Models 

The difficulty in explaining the extended radio halos arises 
from the combination of their ~Mpc size, and the relatively 
short radiative lifetime of the radio emitting electrons. In- 
deed, the difiiusion time necessary for the radio electrons to 
cover such distances is orders of magnitude larger than their 
radiative lifetime. As proposed first by Jaffe (1977), a solu- 
tion to this puzzle would be provided by continuous in situ 
reacceleration of the relativistic electrons. 

An alternative to the reacceleration scenario is given by 
secondary models, which were first put forward by Dennison 
(1980), who suggested that relativistic electrons may be pro- 
duced in situ by inelastic proton-proton collisions through 
production and decay of charged pions. 

There is now increasing evidence that present radio 
data, i.e. the fine radio properties of the radio halos and 
the halo-occurence, may be naturally accounted for by pri- 
mary electron models whereas they may be difficult to be 
reproduced by secondary models (Brunetti, 2003; Kuo et al., 
2004). 

In the framework of primary electron models, merger 



shocks can accelerate relativistic electrons to produce large 
scale synchrotron radio emission (e.g., Rocttigcr ct al., 1999; 
Sarazin 1999; Takizawa & Naito, 2000; Blasi 2001). How- 
ever, the radiative life— time of the emitting electrons difi^us- 
ing away from those shocks is so short that they would just 
be able to produce relics and not Mpc sized spherical radio 
halos (e.g., Miniati et al, 2001). In addition, a number of 
papers (Gabici & Blasi 2003; Berrington & Dermer 2003) 
have recently pointed out that the Mach number of the typ- 
ical shocks produced during major merger events is too low 
to generate non-thermal radiation with the observed fiuxes, 
spectra and statistics. 

a) Re-accclcration of a population of relic electrons by 
turbulence difi^uscd in the cluster volume and powered by 
major mergers is suitable to explain the very large scale of 
the observed radio emission and is also a promising possi- 
bility to account for the fine radio structure of the diffuse 
emission (Schlickciscr ct al.l987; Brunetti et al., 2001a,b). 
In this framework, based on relatively simple assumptions, 
Ohno, Takizawa and Shibata (2002), and Fujita, Takizawa 
and Sarazin (2003) developed specific magneto-turbulent 
models for Alfvenic electron acceleration in galaxy clusters. 
More recently, Brunetti et al.(2004) presented the first self- 
consistent and time-dependent model for the interaction of 
Alfven waves, relativistic electrons, thermal and relativistic 
protons in galaxy clusters. These authors proved that, under 
some physical conditions, radio halos and HXR tails may be 
activated by electron acceleration due to MHD waves in- 
jected during cluster merger events. 

1.4 Why a statistical model ? 

So far two works have investigated the statistics of the for- 
mation of radio halos from a theoretical point of view. 

Enfilin and Rottgering (2002) calculated the radio lumi- 
nosity function of cluster radio halos (RHLF). In a first mod- 
elling, they obtained RHLF by combining the X-ray cluster 
luminosity function with the radio-halo luminosity - X-ray 
luminosity correlation, assuming that a fraction, frh — |, 
of galaxy cluster have radio halos. Then, in a slightly more 
accurate modelling, frh was assumed to be equal to the frac- 
tion of clusters that have recently undergone a strong mass 
increase and the radio halo luminosity of a cluster was as- 
sumed to scale with {l+z)~* (due to the incresing IC-losses). 
In a more recent paper, Kuo et al. (2004) calculated the 
formation rate and the comoving number density of radio 
halos in the hierarchical clustering scheme. The model was 
based on two morphological criteria to define the conditions 
necessary to the formation of radio halos : 1) the cluster 
mass must be greater than or equal to a threshold mass ad- 
justed to observations (Giovannini et al. 1999); 2) the merger 
process must be violent enough to distrupt the cluseter 
core, and thus the relative mass increase was required to 
be Am = {AM /M)th ~ 0.6 according to numerical simu- 
lations (Salvador- Sole et al. 1998). Given the above criteria 
and making use of the Press- Shechter formalism these au- 
thors found that a duration of the radio halo phenomenon 
of the order of 1 Gyr would result to be in good agreement 
with the observed probability of formation of radio halos 
with the mass of the parent clusters. 

As already pointed out, all these approaches are based 
on assumptions in defining the conditions of formations of 
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radio halos based on observational correlations and/or mass 
thresholds. On the other hand, no eflFort has been done so 
far to model the formation of radio halos and HXR tails 
in a self-consistent approach, i.e. an approach which should 
model, at the same time, the evolution of the thermal prop- 
erties of the ICM of the host galgixy clusters and the gener- 
ation and evolution of the non-thermal phenomena. 

As mentioned above, one of the ideas that is produc- 
ing the most promising results for the interpretation of 
non-thermal phenomena in galaxy clusters consists in the 
magneto-turbulent re-acceleration of relic relativistic elec- 
trons leftover of the past activity occurred within the ICM. 

The aim of this paper is thus to obtain a basic modelling 
of the statistical properties of radio halos and HXR tails in 
the above self-consistent approach in the framework of the 
magneto-turbulent re-acceleration scenario. 

In order to have a straightforward comparison with pub- 
lished observational constraints, an Einstein de Sitter (EdS) 
model {Ho — 50 km s^^Mpc"^, qo = 0.5) is assumed in the 
paper. In the Appendix A the results are also compared with 
those obtained in a A-CDM model. 



cluster volume for each merger event, is injected for- and 
thus dissipated in a dynamical crossing time. 
Hi) Particle Acceleration: We focus on the electron compo- 
nent only because the major damping of MS waves (which 
determines the spectrum of these waves and thus the cfh- 
ciency of the particle acceleration. Sect. 5.2) is due to ther- 
mal electrons and thus hadrons cannot significantly affect 
the electron-acceleration process * (see Sect. 5. 2. 2). We as- 
sume a continuous injection of relativistic electrons in the 
ICM due to AGNs and/or Galactic Winds. At each time 
step, given the spectrum of MS waves and the physical con- 
ditions in the ICM, we compute the time evolution of rel- 
ativistic electrons by solving a Fokkcr-Planck eqxiation in- 
cluding the effect of electron acceleration due to the coupling 
between MS waves and particles, and the relevant energy 
losses. 

Given a population of galaxy clusters by combining i)- 
iii) we are thus able to follow in a statistical way the cosmo- 
logical evolution of the spectrum of the relativistic electrons 
in the volume of these clusters and the properties of the 
thermal ICM. 



2 THE STATISTICAL 

MAGNETO-TURBULENT MODEL: 
OUTLINE 

In this Section we outline the formalism and procedures used 
to develope our statistical model. The major steps can be 
sketched as follows : 

i) Cluster formation: The evolution and formation of 
galaxy clusters is computed making use of a semi-analytic 
procedure based on the hierarchical Press & Schechter 1974 
(PS) theory of cluster formation. Given a present day mass 
and temperature of the parent clusters, the cosmological 
evolution (back in time) of the cluster properties (merger 
trees) are obtained making use of Monte Carlo simulations. 
A suitable large number of trees allows us to describe the 
statistical cosmological evolution of galax:y clusters. 

a) Thirbulence in Galaxy Clusters: The turbulence in 
galaxy clusters is supposed to be injected during cluster 
mergers. 

The energetics of the turbulence injected in the ICM is 
calibrated with the PdV work done by the infalling sub- 
clusters in passing through the volume of the most massive 
one; it basically depends on the density of the ICM and 
on the velocity between the two colliding subclusters. The 
sweaped volume in which turbulence is injected is estimated 
from the Ram Pressure Stripping (e.g., Sarazin 2002 and ref. 
therein). We assume that a relatively large fraction of the 
turbulence developed during these mergers is in the form 
of fast magneto-acoustic waves (MS waves). We use these 
waves since their damping rate and time evolution basically 
depend on the properties of the thermal plasma which are 
provided by our merger trees for each simulated cluster. 

The spectrum of the MS waves depends on many unknown 
quantities thus we adopt two extreme scenarios : the first one 
assumes a broad band injection of MHD waves (Sect. 4.2), 
the second one assumes that turbulence is injected at a single 
scale (Appendix B). In both cases the spectrum of MS waves 
is calculated solving a turbulent-diffusion equation in the 
wavenumber assuming that the turbulence, injected in the 



3 CLUSTER FORMATION: PS FORMALISM 
AND A MONTE CARLO APPROACH 

3.1 PS formalism and Merger rate 

The PS theory assumes that galaxy clusters form hierarchi- 
cally via mergers of subclusters which develope by gravita- 
tional instability of initially small amplitude Gaussian den- 
sity fluctuations generated in the early Universe. Making use 
of the PS formalism, it is possible to obtain the probabil- 
ity that a "parent" cluster of mass Mi at a time t\ had a 
progenitor of mass in the range M2 — > M2 + dM2 at some 
earlier time t2, with Mi > M2 and ti > t2. This is given by 
(e.g., Lacey & Cole 1993, Randall, Sarazin & Ricker 2002): 



'P{M2,t2\Mi,ti)dMi 



1 Ml 5c: 



2^M2 (ct2-Ct2)3/2 

{Sc2 — Scl) 



dM2 



exp 



2{ai - ai) 



rfMo 



(1) 



where 6c{z) is the critical linear overdensity for a region to 
collapse at a redshift z; for the EdS model it is given by: 



Sc{t) 



3(127r) 



2/3 



to \ 2/3 



20 



■(7) 



(2) 



with to the present time, cr(M) is the rms density fluctuation 
within a sphere of mean mass M. In Eq. (1) it is 5ci = 5c{t = 
ti) and o"! = a (Ml), with similar definitions for Sc2 and (72. 
The standard deviation of matter density fluctuations at the 
smoothing scale R (cr(M)) is given by (e.g., Peebles, 1980): 



a(M) 



(277)3 



Pf{k)W^{kRM)d'^k 



(3) 



* This is different from the case of Alfven waves whose damping 
may be indeed dominated by the presence of relativistic hadrons 
(Brunetti et al.20O4). 



© 0000 RAS, MNRAS 000, 000-000 



4 R. Cassano & G. Brunetti 



where W^{kRM) is the Fourier transform of the window 
function, the scale Rm is chosen to contain a mass M, and 
P{k) is the matter power spectrum at a given redshift that 
can be expressed as: 



Pf{k) = Pi{k)T^{k;Zf) 



D{zf) 



oc k'^r{k;zf) 



D{zf) 



Dizi] 



(4) 



where n is the spectral index of the primordial power spec- 
trum Pi{k) oc fe"; T{k) is the 'transfer function' which trans- 
fer the initial perturbation to the present epoch and D{z) is 
the growth factor of the linear perturbation. Over the range 
scales of our interest it is sufficient to consider a power- 
law spectrum of the density perturbation given by (Randall, 
Sarazin & Ricker 2002): 



ct(M) = as 



f-V 



(5) 



where ag is the present epoch rms density fluctuation on a 
scale of 8 h'^ Mpc, and Ms = (47r/3)(8ft"^ Mpc)'^p is the 
mass contained in a sphere of radius 8 h^^ Mpc (p is the 
present epoch mean density of the Universe) . The exponent 
in Eq.(5) is given by (Bahcall & Fan 1998) a = {n + 3)/6, 
n being the slope of the power spectrum of the fluctuations. 
Following Randall, Sarazin & Ricker (2002), it is as = 0.514 

for the EdS models. ^ 

It is convenient (Lacey & Cole 1993) to replace the mass 
M and time t (or redshift z) with the suitable variables S = 
a^{M) and x = dc{t). With these definitions, 5* decreases 
as the mass M increases, and x decreases with increasing 
cosmic time t. 

Let K{AS, Ax) dAS be the probability that a cluster 
had a progenitor with a mass corresponding to a change in 
S oi AS = ai - al in the range AS ^ AS + dAS at an 
earlier time corresponding to Ax = 5c2 — <5ci. Then, from 
Eq.(l) one has : 



K{AS, Ax)dAS : 



Ax 



2^ (A5)3/: 



■ exp 



(Aa 



2A5 



3.2 Monte Carlo Technique and Merger Trees 

Following a relatively standard procedure adopted in the 
literature (e.g., Randall et al., 2002; Gabici & Blasi 2003), 
we employ a Monte Carlo technique to construct merger 
trees. Each tree starts at the present time with a cluster of 

mass M and temperature T. We step each simulated cluster 
back in time, using a small but finite time step corresponding 
to a positive increase Aa;. The step size determines the value 
of the minimum mass increment of the cluster, AMc, which 
is due essentially to a single merger event (Lacey & Cole 
1993) : 



{Axf <, 



dlna^ 



rflnM 



/AMc 
V M 



S, 



(7) 



t By comparing Eq. 5 with numerical solutions and analytical 
fits of Eq. 3 (kindly provided by G.Tormen) we find that, for 
the typical masses of the merging subclumps which dominate the 
injection of cluster turbulence in the ICM (> W^^Mq, Sect. 4), 
Eq. 5 is appropriate within a few percent. 
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^0.6 



V 
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0.2 



-5 -4 -3 -2 -1 

log(AS) 

Figure 1. Comulative probability distribution V{< AS, Ax) as 
a function of log{AS). 



where M is the mass of the cluster at the current time 
step. The value AMc gives the mass of the smallest merging 
subcluster we can resolve individually in our trees; we choose 
AMc — W^^h~^ Mq. Thus mass increments smaller than 
this value are considered to be part of the continuous mass 
accretion process in galaxy clusters. 

In order to follow the probability that a merger with a 
given AS (i.e. AM) occurs at a given time we make use of 
the cumulative probability distribution of subcluster masses: 



dAS . (6) V{< AS, Ax) : 



Jo 



K {AS', Ax) dAS' = 



(8) 



whore crfc() is the complementary error function. The 
cumulative probability distribution (Fig. 1) is defined such 
that V{< AS, Ax) 1 for AS ^ oo. 

The Monte Carlo procedure selects a uniformly- 
distributed random number, r, in the range 0-1, then it 
determines the corresponding value of AS solving numeri- 
cally the equation 7'(< AS, Aa;) = r (Fig. 1). The value of 
S2 of the progenitor is given by S2 = Si + AS. The mass of 
one of the subclusters is given by a^{M2) = S2, where a{M2) 
is given by Eq.(5), wheras the mass of the other subcluster 
is AM = Ml - A/2. 

We define Mmin = min(M2,AM) and Mmax = 
max(M2, AM). In order to speed up the computational pro- 
cedures, without significantly affecting the results, we con- 
sider two cases : 

• i) if Mmin < 1 X 10^'^ Mq and Mmin < AMc, then the 
event cannot be identified with a merging subcluster and is 
considered as accretion. If Mmin < 1 x 10^^ Mq and Mmin > 
AMc, then the event is considered a very minor merger and 
its contribution to the injection of cluster turbulence (Sec. 4) 
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is neglected ^ . In both cases, the mass of the parent cluster is 
simply reduced to M2 = Mi — Mmin and the next time-step 
in the merger tree starts from M2. 

• ii) if Mmin > 1 X 10^^ Mq then the event is treated as 
a merger and we calculate all the physical quantities useful 
for the computation of the energy of the turbulence gener- 
ated during this event (Sect. 4). In this case, if Mmin is also 
greater than a given value of interest (Sect. 7) we follow back 
in time the evolution of both the subclusters (i.e., Mmin and 
Mmax) constructing the merger tree for each subcluster. 

This procedure is thus iterated until either the mass of 
the larger cluster drops below AMc or a maximum redshift 
of interest Zmax is reached. An example of a merger tree 
obtained from our procedure (tracing the evolution of the 
Mmax clusters only) is shown in Fig. 2 as a function of both 
look back time and redshift. 

Our procedure is basically a Binary Merger Tree Method 
which does not allow to describe multiple nearly simultane- 
ous mergers. This simple procedure, however, is sufficient 
for our purposes since multiple mergers mainly afi'ect the 
evolution of low mass halos at relatively high redshift which 
are not interesting for the study of the non-thermal phe- 
nomena. The implementation of more complicated N-Branch 
Tree Methods can be found in Somerville & KoUatt (1999). 



4 RAM PRESSURE STRIPPING, 

TURBULENCE AND MHD WAVES 

4.1 Turbulence injection rate 

The passage of the infalling subhalos through the main clus- 
ter during mergers induce large-scale bulk flows with veloc- 
ities of the order of ~ 1000 km s~^ or larger. Numerical sim- 
ulations of merging clusters (e.g., Roettiger, Loken, Burns 
1997; Ricker & Sarazin 2001; Tormen et al. 2004) provide 
a detailed description of the gasdynamics during a merger 
event. It has been found that subclusters generate laminar 
bulk flows through the sweeped volume of the main clus- 
ters which inject eddies via Kelvin-Helmholtz instabilities at 
the interface of the bulk flows and the primary cluster gas. 
Finally these eddies redistribute the energy of the merger 
through the cluster volume in a few Gyrs by injecting ran- 
dom and turbulent velocity fields. 

The impact velocity between the subclusters increases 
at the beginning of the merger and then it saturates when 
the subclusters interpenetrate each other. Depending on the 
initial conditions and on the mass ratio of the two subclus- 
ters, during the merging process the infalling halos may be 
efficiently stripped due to the ram-pressure. However, the 
numerical simulations show that the efficiency of the ram 
pressure stripping is reduced by the formation of a bow 
shock on the leading age of the subcluster. This bow shock 
forms an oblique boundary layer which slows the gas flow 
and redirect it around the core of the subcluster so that, 
at least in the case of mergers with mass ratios < 10, a 
significant amount of the subcluster gas is found to be still 



"f- Note that we are interested in describing mergers of typically 
> 5 X IO^^Mq, Sect. 7. 
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Figure 2. Example of Merger Trees obtained from Monte Carlo 
simulation in a EdS universe for clusters with present day mass 
Afo = 2.5 X W^Mq: a) Log{M) - z; b) M/Mq - t/*0 with t 
cosmic time; to present time. 



self-bounded after the first passage through the central re- 
gions of the main cluster (Roettiger et al. 1997; Tormen et 
al. 2004). 

Due to the complicated physics involved in these events, 
the details of the injection and evolution of turbulent mo- 
tions in galaxy clusters during merging processes are essen- 
tially a still unexplored issue. However, turbulence should be 
basically driven by the PdV work done by the infalling halos 
through the volume of the primary cluster and the turbu- 
lent motions should be initially injected within the volume 
sweeped by the passage of the subhalos (e.g., Fujita, Tak- 
izawa, Sarazin 2003). Following this simple scenario, in this 
Section we estimate the rate of turbulence injected during a 
merger event. As a necessary approximation (due to the PS 
formalism) in the calculations, we assume that subclusters 
undergo only central collisions. 

The relative impact velocity of two subclusters with 
mass Mmax and Mmin which collide (at a distance Rmax be- 
tween the centers) starting from an initial distance do with 
zero velocity is given by (e.g., Sarazin 2002): 




(9) 
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Figure 3. Panel a): Ratio between the energy injected in form of turbulence and the thermal energy of the system as a function of 
the mass ratio of the two subclusters. Thermal energy is calculated for both the sum of the two subclusters (lower lines) and for the 
main cluster alone (upper lines). We stress that in the case of a mass ratio < 1.5 our approach is quite inadequate because essentially no 
primary cluster exists, however these mergers are very rare events and do not dominate the injection of turbulence in our model. 
Panel b); Ratio between the stripping radius and the radius of the radio halos (assuming Rjj = 500h~^ kpc) reported as a function of 
the mass ratio of the two subclusters. 

In both panels calculations are obtained following the recipes given in Sect. 4.1 for a M = 5, (solid lines), 2 (dotted lines), and 1 X IO^^Mq 
(dashed lines) clusters. 



where do — rjvRn 



Rmax is the virial radius of the main cluster, i.e. 

1/3 



Rmax 



4n Ac{z)pcrit{z) 



and 



(10) 



pcrit{z) is the critical density at redshift z and Ac{z) = 
IStt^ ~ 180 is the ratio of the average density of the cluster 
to the critical density at redshift z. 

While the smaller subcluster crosses the larger one, it is 
stripped due to the effect of the ram-pressure. The strip- 
ping is efficient outside a radius (stripping radius) at 
which equipartition between static and ram-pressure is es- 
tabhshed, i.e. : 



2 Pmin 



/imp 



(11) 



where, as an approximation, p^^^ is fixed at the average 
density of the ICM of the larger subcluster : 



Pmax 



ax 



X fb 



(12) 



with ft = 0.25 (^) ^^'^ the observed barion fraction of clus- 
ters (Ettori & Fabian 1999, Arnaud & Evrard 1999). Eq.(ll) 
is solved numerically at each merger event assuming that the 
density profile of the ICM of the smaller cluster, pmin, is de- 
scribed by a /3-model (Cavaliere & Fusco-Femiano, 1976): 



pminij'^ — Pmm (0) 



l + (-)^ 



-3/3x/2 



where the normalization is given by 



(0) 



in 



drr 



l+(f)^ 

' C 



-, -3/3x/2" 



(13) 



(14) 



where a core radius Tc = 0.1 Rmin and /3x — 0.8 are as- 
sumed. The temperature of the smaller cluster, Tmm, in 
Eq.(ll) is estimated making use of the observed M-T rela- 
tionship (e.g., Nevalainen et al. 2000). As a general remark 
we stress that the value of the stripping radius obtained 
above would give the mean value of during a merger and 
it is not the minimum r^. In qualitatively agreement with 
numerical simulations, this approach yields rs —> in the 
case of mergers with large mass ratios between the two col- 
liding subclusters. 

The motion of the smaller cluster through the ICM of the 
main one generates fiuid turbulence. Following Fujita et 
al.(2003) we assume that turbulence is initially injected in 
the sweeped volume, Vt ~ TvrlRmax, with a maximum tur- 
bulence lenght scale of the order of 2 x r^. The total en- 
ergy injected in turbulence during a merger event is thus 
■Et Pmax,s'"?Vt, where p„^^_^ is the ICM density of the 
main cluster averaged on the sweeped cylinder. We assume 
that the duration of the injection is of the order of a crossing 
time, Tcros — Rmax/vi, then the turbulence is dissipated in 
a relatively short time. 

The use of the averaged density of the ICM of the primary 



© 0000 RAS, MNRAS 000, 000-000 



Cluster mergers and non-thermal phenomena : a statistical magneto-turbulent model 7 



64 



^ 63 



O 



62 - 



61 



I I i r- 

M =lxlOi5M^ 




0.2 0.4 
z 



0.6 



O.B 



Figure 4. Evolution of the thermal energy (solid lines) and of 
the energy injected in fluid turbulence (dashed lines) integrated 
during cluster life (at redshift z) in typical galaxy clusters. The 
thin lines are for a cluster with present time Mq = 10^^ Mq and 
the thick lines are for a cluster with present time Mo = 5 x 



cluster, p^axj of the initial impact velocity between the sub- 
clusters, Vi, and of the density of the main cluster averaged 
on the sweepod cylinder, pmaa: s ^^'^ caluclations of the in- 
jected turbulence is a necessary simplification which however 
guarantees a basic estimate of the averaged injected turbu- 
lence in the ICM and which does not depend on essentially 
unknown details. For seek of completeness, in Fig. 3a we 
report the tj'pical ratio between turbulent energy injected 
by a merger event and the thermal energy of the system as 
a function of the mass ratio between the two colliding sub- 
clusters; it is found that major mergers may channel about 
10-15 % of the thermal energy in the form of large scale tur- 
bulence. In Fig. 3b we also report the value of the stripping 
radius as a function of the mass ratio of the two colliding 
subclusters. It is found that Vs (i.e., the mean value of 
during a merger event) is typically larger than the radius of 
the radio halos, Rh, for the merger events which mainly con- 
tribute to the injection of cluster turbulence in our model. 
If the sweeped volume is smaller than that of the radio halo, 
we assume that the injected turbulence is diffused over the 
volume of the radio halo, Vh = ^ttRh, which is basically 
equivalent to assume that the integral cross section of the 
ensemble of minor mergers which occur in a time interval of 
~Gyr is comparable to Rh- 

Under these hypothesis, the injection rate per unit vol- 
ume of turbulence is given by : 



Et 



X Vh 



Pmax.s 3 
— —Vi 



Vh 



(15) 



As a relevant example, in Fig. (4) we report the cos- 
mological evolution of the thermal energy of galaxy clusters 
with different masses, together with the total energy injected 



in form of turbulence in the ICM. The energy in turbulence 
is calculated by integrating the contributions from all the 
merger events. The thermal energy of the considered clus- 
ters, calculated assuming the observed M-T relation (e.g., 
Nevalainen et al. 2000) , increases from about lO'^^erg at 
a ~ 1 to a few lO^'^erg at the present epoch depending 
on the mass of the cluster. As it should be, we note that 
the energy budget injected in turbulence during cluster for- 
mation is well below the thermal energy; this indicates the 
consistency of our calculations. In particular the turbulent 
energy is found to be ~ 15% that of the thermal energy in 
agreement with recent numerical simulations (Sunyaev et al. 
2003) and with very recent observational claims (Schuccker 
et al. 2004). Finally, as reasonably expected, the energy in- 
jected in turbulence calculated with our approach is found 
to roughly scale with the thermal energy of the clusters. 



4.2 Spectrum of the magnetosonic waves 

Cluster mergers should generate magnetosonic (MS) waves 
in the ICM, the injected energy and spectrum of these waves 
depends on many unknown quantities. A reasonable attempt 
is to assume that a fraction, r]t, of the energy of the turbu- 
lence (Sect. 4.1) is in the form of MS waves. We shall consider 
two extreme scenarios: 

i) in the first one we assume that MS waves arc driven by 
the plasma instabilities (e.g., Eilek 1979, and ref. therein) 
which develope in the turbulent field generated during clus- 
ter mergers. In this case MS waves may be injected over a 
broad range of scales. Here, we shall adopt a simple power 
law injection spectrum of these waves: I{k) = lok'" for 

k > kmin ^ Tt/^s; 

ii) in the second one, we assume that MS waves are basi- 
cally injected at a single scale, k = kmin ~ T^/rs, from which 
a MHD turbulence cascade is originated. 

In both cases the decay time of the MHD turbulence at 
the maximum/injection scale, Linj ~ 2rs, can be estimated 
as (e.g.. Appendix B) Tkk{Li„j) ~ one has : 



rfcfc(Gyr) ~ 1 x ( 



2 • lO^km/s 



)(^)-^ 
500kpc^^0.25^ 



(16) 



which is of the order of a crossing time and thus allows the 
MHD turbulence to diffuse filling a volume of the order of 
that of radio halos (or larger) with a fairly uniform intensity. 

In the following we focus on the first scenario, while in 
the Appendix B wo consider the second picture. Appendix B 
demonstrates that these two extreme scenarios lead to very 
similar results and thus that, in our model, the details of the 
injection process of the MS waves do not appreciably change 
the conclusions. 

In the case in which a power law spectrum of MS waves 
is injected in the ICM, one has : 



lok °'dk = rjt - 



Et 



(17) 



with kmin ~ Tr/rs and kmax = i^p/vM i^p being the pro- 
ton cyclotron frequency and vm the magnetosonic velocity 
(Eq. 34). From Eq. (17) we find: 
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lo 



VjLf min 



(<^-i) 
v„ 



k: 



ij-s < Rh ) 



(18) 



Thus, the injection of the MS waves is obtained by com- 
bining Eqs.(15) and (18). 774 is the first free parameter of our 
model, in order to have a self-consistent modelling it should 
be rit < 1. 

In general, the spectrum of MHD waves injected in the 
ICM evolves due to wave-wave and wave-particle coupling. 
The combination of these processes produces a modified 
spectrum of the waves, VKk(t). In the quasi linear regime 
the spectrum of the waves can be calculated solving a con- 
tinuity equation in the wavenumber space: 



at 



dk 



' i=l 

+Ht) 



(19) 



The first term on the right hand describes the wave-wave 
interaction, with diffusion coefficient Dkk = k^ /ts (with 
Ts the spectral energy transfer time). The second term in 
Eq.(19) describes the damping with the relativistic and ther- 
mal particles in the ICM, and the damping due to thermal 
viscosity. In the following, we shall neglect the term due 
to the wave wave interaction, this is justified provided that 
the time- scale of the dampings are smaller than that of the 
wave-wave cascade (or comparable), at least for the range 
of scales which contribute to the acceleration process (as- 
suming spectra fiatter than Wk oc k~^; see also Appendix 
B). Under physical conditions typical of the ICM the most 
important damping in the coUisionless regime is that with 
the thermal electrons (e.g., Eilek 1979; see also Sect. 5.2.2). 
An estimate of this damping rate (for va < vm as in the 
ICM) is given by Eilek 1979; a relatively simple formula, 
consistent whitin a 10% with the Eilek's results, is: 



rth,e = V32^ nth (meKBT) ( ^ ^ j^^) fc 



(20) 



where nth is the number density of the thermal electrons, 
is the turbulent magnetic energy density, = Wk x 
(1 -I- 87rP/B^)"\ with P ~ 2nthKBT the thermal pressure, 
B the plasma magnetic field (e.g., Barnes & Scargle 1973) 
and I{x) is a numerical value given by : 



I{x) 



dx 



Va; x'^ J 



(21) 



where x = 



— , with p|| the component of the mo- 
mentum of the thermal electrons along the magnetic field 
lines and vth = {2 Kb T/rrieY^'^. 

Since for each merger event we are interested in the evo- 
lution of the spectrum of the injected waves on a time scale 
of ~ 1 Gyr, which is orders of magnitude longer than the 
typical time scales of the damping processes, the spectrum 
of the waves is expected to approach a stationary solution 
{dW]^/dt = 0). Prom Eq. (19) this solution is given by : 



/(k) 



/(k) 



rth,e(fc) /(T) k 



(22) 



In Sect. 5.2 we will derive the efficiency of electron acceler- 
ation due to the MS waves. Here we would just point out 
that the acceleration time, race, depends on : 

T„„,: oc / kWkdk (23) 



.1 

CO oc / 



which leads (making use of Eqs. 18 and 22) to the nice 
result that the acceleration time in our model, and under 
our assumptions, does not depend on the slope of the in- 
jection spectrum of MS waves (which depends on basically- 
unknown details of the injection mechanism) and on the 

value of kmin- 



4.3 Spectrum of MS waves during cluster 
formation 

In this Section we estimate the spectrum of MS waves re- 
sulting from the combination of the contributions of several 
mergers during the process of cluster formation. For a given 
galaxy cluster, we define to be the redshift at which the 
j"' merger event starts. For an Einstein-De Sitter model, 
the corresponding time, , is : 

= _^ 1 (24) 

' 3ffo (1 -1-^1)3/2 ^ 

In our simple modelling we assume that the duration time 
of a merger is of the order of a crossing time, At = tiros, 
and that turbulence is injected during this time interval and 
then suddenly dissipated via damping processes § . Thus the 
turbulence injected during the j*'' merger is dissipated at 
time fj: =tl + tiross and the corresponding redshift is given 
by: 



4 = i 



3Hot} 



^2/3 _ ^ 



(25) 



We describe the spectrum of the MS-waves established dur- 

■th 



ing the j merger event as : 

Wi{z) = Wi{zi)xS'iz) 

were {z) is a step function defined as : 



1 iz}<z< zf) 
(otherwise) 



(26) 



(27) 



According to the hierarchical scenario adopted in this 
paper, clusters undergo several merger events which con- 
tribute to the injection of turbulence yielding a combined 
spectrum of MS-waves. Since under stationary conditions 
and neglecting the wave wave interaction term, Eq.(19) is a 
linear difi'erential equation, the spectrum of the MS-waves 
resulting from the combination of the different merger events 
is given by the sum of all the contributions (Eq.26), i.e. : 



(28) 



§ Note that if the injection time is slightly longer than tcros then 
the probability to combine the effect of several mergers increases 
and the efficiency of the model would slightly increase (see above, 
Eq.28). 
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5 PARTICLE EVOLUTION AND 
ACCELERATION 

Volk et aJ. (1996) and Berezinsky et al. (1997) have shown 
that cosmic ray protons are most likely accumulated in 
galaxy clusters as their diffusion time scale is of the order or 

larger than the Hubble time. An oven stronger conclusion 
can be applied to the case of the relativistic electrons in- 
jected in the ICM. Indeed, by assuming a Kolmogorov spec- 
trum for the turbulent field the parallel diffusion length, Lo, 
of relativistic electrons is given by (e.g., Brunetti 2003): 



120 



7/10' 

B^G 



lo 



TD 



100 kpc I Gyr 



(kpc) 



(29) 



whore lo is the maximum coherence scale of the mag- 
netic field fluctuations. We notice that 7 ~ 100 — 300 relic 
relativistic electrons cannot diffuse more than 200 kpc dur- 
ing their life-time (a few Gyrs, e.g., Sarazin 1999). When 
turbulence is generated in galaxy clusters the ensuing in- 
crease of the particle-wave scattering frequency would fur- 
ther reduce the net diffusion coefficient. Only in the case of 
very strong turbulence anomalous particle turbulent diffu- 
sion may allow electrons to diffuse over larger scales (e.g., 
Duffy 1994; Enfilin 2003), but the requested conditions are 
well beyond those assumed in the present paper. 
Thus, we can safely assume that electrons injected by some 
mechanism in the ICM simply follow the thermal plasma 
and magnetic field. Under this condition, the time evolution 
of relativistic electrons with isotropic momentum distribu- 
tion is provided by a Fokker-Planck equation (e.g., Tsytovich 
1966; Borovsky & Eilek 1986) for the electron number den- 
sity : 



dN{p,t) 
dt 



dp 



N{p,t) 



dp I 

dt 'rad 
d 



dp 



I '^P I 

dN{p,t) 



dp 



+ Qe{p,t) (30) 



where Dpp is the electron diffusion coefficient in the 
momentum space due to the interaction with the MS waves, 
dp/dt^ and dp/dt^^ are the terms due to ionization and ra- 
diative losses and Qe is an isotropic electron source term. 



5.1 Particle injection 

A relevant contribution to the injection of cosmic rays 
in clusters of galaxies comes from Active Galactic Nuclei 
(AGN). AGNs indeed inject in the ICM a considerable 
amount of energy in relativistic particles and also in mag- 
netic fields, likely extracted from the accretion power of their 
central black hole (Enfilin et al., 1997). 

Powerful Galactic Winds (GW) can inject relativistic 
particles and magnetic fields in the ICM (Volk & Atoyan 
1999). Although the present day level of starburst activity 
is low and thus this mechanism is not expected to produce 
a significant contribution, it is expected that these winds 
were more powerful during early starburst activity. Some 
evidence that powerful GW were more frequent in the past 
comes from the observed iron abundance in galgixy clusters 
(Volk et al. 1996). 

In addition, cluster formation is also believed to provide 



a possible contribution to the injection of cosmic rays in the 
ICM due to the formation of shocks which may accelerate 
relativistic particles (Blasi 2001; Takizawa & Naito, 2000; 
Miniati et al., 2001; Fujita & Sarazin 2001). The efficiency 
of this mechanism is related to the Mach number of these 
shocks which is an issue still under debate. Semi-analytical 
calculations based on PS-Monte Carlo techniques (Gabici & 
Blasi 2003; Berrington & Dermer 2003) find that the bulk of 
the shocks have Mach numbers of order 1.5, as also observed 
by Chandra (e.g.. Markevitch et al. 2003). First cosmolog- 
ical numerical simulations found that the Mach number of 
merger and internal flow shocks peacks at ~ 5 (Miniati 
et al. 2000), and that the bulk of the energy is dissipated 
at 4 < < 10 shocks (Miniati 2002). On the other hand, 
more recent numerical simulations find more wealc shocks 
with the bulk of the energy dissipation (thermal energy and 
cosmic rays) at internal shocks with 2 < M < A (Ryu et al., 
2003). However, the comparison with semi-analytical calcu- 
lations appears difficult because of a different classification 
of the shocks in the two approaches. 

Independently from the specific scenario adopted for the 
injection of relativistic particles, a power law spectrum for 
the injection rate of relativistic electrons up to a maximum 
momentum, Pmax, can be reasonably assumed : 



Q,{p,t) = K,{t)p- 



(31) 



We parameterize the injection rate by assuming that the 
total energy injected in cosmic ray electrons (for p > pmin) 
during the cluster life is a fraction, rje, of the total thermal 
energy of the cluster at ^ = 0, i.e. : 

ne = — dr Qeip,T)pdp (32) 

where fth is the present day thermal energy density of 
the ICM. 

The injection rate should depend on the number and en- 
ergetics of AGNs and GWs in galaxy clusters which are ex- 
pected to be considerably larger at high redshifts. However, 
since electrons injected at relatively high redshifts cool very 
rapidly because of the combination of high energy losses and 
low efficiency of the particle acceleration mechsnism (Fig. 5), 
only the electrons injected at relatively low redshifts can 
be accelerated and therefore contribute to the non-thermal 
emission observed at z < 0.2. As a simplification, we adopt 
a constant injection rate of electrons so that (for s > 2) the 
normalization of the spectrum of the injection rate is given 
by: 

(33) 



where th is the Hubble time, rje is the second free parameter 
in our model. In the following we use s = 2.5, Pmin/mc = 60, 
and pmetx/mc = 10*; as discussed in Sect. 8 the basic results 
of our model do not depend on the values adopted for these 
parameters. 



5.2 Energy gains and losses of the relativistic 

particles 

Fast MS waves can accelerate relativistic particles via 
particle-wave interaction. This interaction occurs with both 
thermal (expecially) and relativistic particles via Landau 
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damping, the strongest damping of MS waves in the col- 
lisionless regime. The necessary condition (Melrose 1968; 
Eilek 1979) is a; — = 0, where uj is the frequency of 
the wave, is the wavenumber projected along the mag- 
netic field, and — v/i is the projected-particlc velocity. 
Under the typical conditions of the ICM, the dispersion re- 
lation for MS waves in an isotropic plasma is given by (e.g., 
Eilek 1979): 



■ Vm ^ -V, 



2,2 
ion + f'A 



(34) 



fe2 "™ - 3 ■ 

where Vion and va are the ion-sound velocity and the Alfven 
velocity, respectively. In a isotropic distribution of waves and 
particle momenta the diffusion coefficient in the momentum 
space, for va < vm, is given by Eilek (1979). This can be 
expressed as: 



I)pp(p,t)~4.45 



2 2 pkmax 
2Vm P / 



kWk{t)dk 



(35) 



The acceleration time scale, which in this case does not 
depend of the particle energy, is given by : 



'''acc — X — 2 



pp 

,2 



(36) 



and thus the systematic energy gain of particles interacting 
with MS waves is given by : 

ays 

■ XP (37) 



\dtj 

acc 



5. 2. 1 Electrons 

Relativistic electrons with momentum pe = ruecy in the 
ICM lose energy through ionization losses and Coulomb col- 
lisions (Sarazin 1999): 



= -3.3 X 10" 



nth 



1 + 



ln(7/wth) 
75 



(38) 



where nth is the number density of the thermal plasma. 

Relativistic electrons also lose energy via synchrotron 
emission and inverse Compton scattering off the CMB pho- 
tons: 



it) 



rad 



-4.8 X 10 



-4 2 
P 



my 



2/3 



(39) 



larger than ~ 300 MeV. The time scale of this process is 
Tpp ~ 30(nth/10"^)"^Gyr. The interactions with the ICM 
may generate an appreciable flux of gamma rays and neu- 
trinos, in addition to a population of secondary electrons 

(Blasi & Colafranccsco 1999). 

Protons which are more energetic than the thermal elec- 
trons lose energy due to Coulomb interactions. For trans- 
relativistic and sub-relativistic protons this channel can eas- 
ily become the main channel of energy losses in the ICM. 

In our calculations wc focus on the population of rela- 
tivistic electrons, do not consider proton acceleration, and 
neglect the effect of these particles on the efficiency of the 
electron acceleration. 

The resonance condition VMk = fc||t;|| implies that only 
a very small fraction of MS waves (those making an angle 
~89-91 degrees with the local B-field) cannot be damped 
by the thermal electrons, but only by the relativistic parti- 
cles (protons and electrons), while outside this narrow cone 
the damping due to the thermal electrons should be the 
strongest one (e.g., Eilek 1979). As a consequence, since 
in our calculations we assume a continuous pitch angle 
isotropization (e.g.. Miller et al. 1996) and an isotropic dis- 
tribution of MS waves which propagate in a complex geom- 
etry of the field lines, the damping of MS waves should be 
dominated by the effect due to the thermal electrons in the 
ICM. 

As an example, we assume that MS waves are injected in 
the central -^l Mpc'^ of a massive cluster for 0.5 T Gyr with 
a total energy budget of the order of that of the thermal 
ICM within the same region. We calculate particle acceler- 
ation and find that about ~ 4 — 10% of the energy flux of 
these waves is channelled into the acceleration of relativistic 
protons (assuming an initial energy density of these parti- 
cles of the order of few % of the thermal energy density and 
Np{p) oc p"^'^): this corresponds to < 1% of the total ther- 
mal energy of the cluster. At the same time, wo find that 
only ~ 0.1% of the energy flux of the MS waves should be 
channelled into the acceleration of relativistic electrons to 
produce an HXR luminosity of ~ 10"'^ erg s~^ from the same 
volume (r?e ~ 0.003, Fig. 6). Consequently the 90 - 95% of 
the energy flux of the MS waves is channelled into the ther- 
mal electrons and thus the resulting spectrum of these waves 
may be estimated with good approximation by Eq.(22). 

Detailed time dependent calculations which include 
electron and proton acceleration due to MS waves and a 
comparison with the case of Alfven waves will be given in a 
forthcoming paper (Brunetti et al., in prep.). 



where i3^G is the magnetic field strength in nG and is 
the pitch angle of the emitting electrons; in case of efficient 
isotropization of the electron momenta it is possible to aver- 
age over all possible pitch angles, so that < sin^ 9 >= 2/3. 

Eqs.(35), (38) and (39) gives the coefficient of the 
Fokker-Plank equation (30). 

5.2.2 Protons 

For seek of completeness in this Section we briefly discuss 

the case of relativistic protons. 

The main channel of energy losses for these particles is rep- 
resented by inelastic proton-proton collisions, which is a 
threshold reaction that requires protons with kinetic energy 



6 RADIO HALOS AND HXR TAILS 

6.1 Cluster evolution and electron spectrum 

In this Section we combine the formalism developed for the 
evolution of the turbulence (Sects. 3 & 4) with the recipes 
for particle acceleration and evolution (Sect. 5) to model 
the cosmological evolution of the spectrum of the relativistic 
electrons in galaxy clusters. 

The electron-acceleration coefficient, due to the effect 
of MS waves at redshift z, is obtained by combining Eq. (36) 
with Eqs. (35, 18, 22, 28): 
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Figure 5. Top panel: evolution with redshift of the electron- 
acceleration coefficient due to MS waves as obtained from Eq.40 
for a cluster of Mq = 1 x 10^^ Mq at the present time. Bot- 
tom panel: electrons spectra (in arbitrary units) calculated at 
different redshifts (also marked in the top panel) for the same 
cluster. Calculations are performed for s = 2.5, B = 0.5/iG, and 
77i = 0.26. 



2.23 X 10~^'^r?t 
(i?H/500kpc)3 

(r,/500kpc)^ 
(fcT/7keV)i/2 



E 



M^a^ + Mmin 2.6MpC \3/2 
2 X 1015 Mq i?™, J 

1 if < Rh ^ 

{RH/r-sf iirs>RH 



(40) 



where only mergers which contribute to the turbulence 
spectrum at redshift z (Sect. 4.3, Eq.28) are considered. The 
evolution of the electron spectrum is thus obtained from the 
numerical solution of the Fokker-Planck equation (Eq. 30) 
by adopting the values of the coefficient Dpp (Eq. 36 and 
Eq. 40) and of the energy loss terms (Eqs. 38-39) at each 
redshift. 

In Fig. (5) we report an example of the time evolution 
of the electron-acceleration coefficient obtained for a typical 
massive cluster (top panel) and the corresponding spectra of 
the electrons at different relevant times (bottom panel): an 
increase of the acceleration coefficient produces an increase 



of the maximum energy of the electrons. The reported re- 
sults indicate that cluster-merger activity at low redshift 
can generate an increse of the cluster turbulence which may 
be sufficient to accelerate electrons up to 7 >> 10^, neces- 
sary to produce synchrotron radiation in the radio band. 
It should be noticed that electrons are accelerated (and 
cool) with a delay time (of the order of the correspond- 
ing electron-acceleration time ~ X~^) with respect to the 
abrupt increases (decreases) of the values of the accelera- 
tion coefficient. 

In Fig. (6) we show the broad band non-thermal emis- 
sion (synchrotron and IC) from the galaxy cluster reported 
in Fig. 5 assuming rje = 0.003 and B — 0.5/i G (for the sake 
of completeness synchrotron and IC equations are given in 
Appendix C). The aim of this Figure is to show that syn- 
chrotron (and IC) luminosities of the order of those of the 
most luminous radio halos can be reasonably obtained. On 
the other hand, it should be stressed that the synchrotron 
spectrum reported in Fig. (6) is obtained assuming a con- 
stant value of B through the cluster volume. More reason- 
able calculations should assume a radial gradient of the mag- 
netic field strength which causes a stretching in frequencies 
of the synchrotron spectral shape with respect to that of 
Fig.(6) (e.g., Brunetti et al. 2001; Kuo et al. 2003). 

The synchrotron emitted power from radio halos is also 
expected to increase with increasing the mass of the parent 
clusters. Indeed the bolometric synchrotron power roughly 
scales as Psyn oc B^'y^NeRij where 7t, = ^ is the maximum 
energy of the accelerated electrons (/3 is the total energy-loss 
coefficient, Eq. 39) and A'c is the number of relativistic elec- 
trons in the cluster emitting volume. During major mergers, 
from Eq.(40, with T oc M"", a ~ 0.55 - 0.67) one has : 



PsynxN,M''-''g{rs,RHy 



B^ 



(41) 



where g{rs,RH) {g = for < Rh and g — 
rs > Rh) is a slightly increasing function of cluster mass. 
Thus assuming B < 3/iG one can find that the synchrotron 



power would roughly scale as Psyn oc M 



b(M) 



< NcB^ >(M), 



where b{M) ~ 1.5 - 2.5 (& = 1.5 for M > 3 ■ IO'-^Mq and 
b = 2.5 for M < IO^^Mq), and < N^B^ >(m) is expected 
to increase with cluster mass. Although studies of the syn- 
chrotron power - mass (or temperature) correlation for ra- 
dio halos use the monochromatic synchrotron emission at 
1.4 GHz (e.g., Liang et al. 2000; Govoni et al. 2001; Giovan- 
nini & Feretti 2002), the above expectations seem to be in 
line with observations. 



6.2 On the required values of % and r/^ 

In this Section we derive the range of values of the two free 
parameters of our model, r/t and rje, which provide a reason- 
able agreement with the general properties of radio halos. 
In order to check the reliability of the obtained values, these 
are then compared with independent findings and general 
expectations from both analytical and numerical calcula- 
tions. 

The first free parameter is rjt which is defined as the 
fraction of the fluid turbulence in MS waves. The value of 
rjt drives the efficiency of the electron acceleration and thus 
the resulting maximum energy of electrons, 76 = 4, (/3 is 
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Figure 6. Broad band synchrotron (SYN) and Inverse Compton 
(IC) spectra calculated for the case reported in Fig. 5 and for 
Tje = 0.003, Rh = 500 kpc, and B = 0.5^G. Horizontal bars give 
the radio (used for radio halos) and HXR observational range of 
frequencies. 



the total energy-loss coefficient, Eq. 39), and the maximum 
synchrotron emitted frequency fb = ^;|^76- Under the 
assumption that the losses of the electrons are dominated 
by the IC mechanism, the acceleration coefficient is thus 
related to the break frequency by : 



X ■ 



: 6.3 X 10 B„r/ (1 + z) 



(42) 



The values of rjt are constrained by requiring that the 
accelerated electrons can produce synchrotron radiation in 
the radio band with the spectral shape observed in the case 
of radio halos, i.e. with spectral index a = 1.1 — 1.5 be- 
tween 327 and 1400 MHz (e.g., Kempner & Sarazin 2001). 
The synchrotron spectral index between two fixed frequen- 
cies depends on the value of vt and also on the shape of the 
spectrum of the emitting electrons. Given the typical shape 
of the spectrum of the emitting electrons accelerated during 
cluster mergers in our calculations, we are able to estimate 
the minimum typical value of Vb necessary to account for the 
spectral indices of the observed radio halos: Ub > 70 MHz is 
obtained. From Eq. (42), this limit translates into a limit on 
X (given in Eq.40) : 



> Xm«n = 7.4 X 10-^^(1 



0.5 



-1/2 



(43) 



Radio halos have a typical size Rh ~500 Kpc and they are 
found in massive galaxy clusters (M>1O^'^M0). Thus we de- 
rive the value of x for these typical clusters in our synthetic 
population and find that rjt = 0.2 — 0.3 is required to satisfy 
the condition of Eq. (43) during major mergers (at z < 0.2; 
B ^ 0.5/iG is adopted). This is the first important result of 
our modelling since it basically proves that if a fraction of 
the kinetic energy of cluster mergers is channelled into MS 
waves then this is sufficient to power particle acceleration in 



the ICM with the efficiency requested in the case of radio 
halos. Although there are no numerical studies which are 
aimed at a detailed investigation of the cluster turbulence 
injected during merging processes, a general finding of high 
resolution numerical simulations is that a relevant fraction 
(10-30 %) of the thermal energy in galaxy clusters is in the 
form of compressible plasma turbulence (e.g., Sunyaev et al. 
2003 and ref. therein). This is in line with the requirements 
of our model. 

The second free parameter in our model is rj^ which 
gives the ratio between the energy injected in relativis- 
tic electrons during the cluster life and the present day 
thermal energy of the ICM. The values of i]e can be con- 
strained by requiring that the model reproduces the typical 
radio (Lr) and hard-X ray (Lhx) luminosities observed in 
galaxy clusters: Lr = 10''° - 10**'ergs"^ (Feretti 2003) and 
Lhx = 10"*^ - lO^-^ergfi-i (Fusco-Femiano et al.2003). We 
derive the requested values for typical massive galaxy clus- 
ters in our synthetic population during the time intervals in 
which the condition of Eq. (43) is satisfied. Making use of 
Eqs. (33), (CI), and (C3) we find that 7?^ = 10"" - 10"^ 
is sufficient to match the observed luminosities at z < 0.2 
{B ~ 0.5/iG is assumed). The above r^e-values are very rea- 
sonable for massive clusters (e.g.. Sect. 5) and they are also 
much smaller than those assumed in other modellings of 
non-thermal emission from galaxy clusters (e.g., rje ~ 0.1, 
Sarazin 1999). This is mainly because in our model the re- 
sulting spectrum of the emitting electrons during an effi- 
cient acceleration period is not a simple power law, but it 
is peaked at the energies required to emit the synchrotron 
and IC radiation (e.g.. Fig. 6) and this strongly increases 
the emitting efficiency (see also Sect. 8). 



7 STATISTICS AND COMPARISON WITH 
OBSERVATIONS 

In Sect. 6.2 we have derived a criterion for radio halo for- 
mation: clusters may have radio halos if xiVt) > Xmin- By 
making use of this criterion, the goal of this Section is to 
calculate the formation probability of radio halos with clus- 
ter mass and to compare expectations with observational 
constraints. 

In order to have a prompt comparison with observations 
we calculate the formation probability in the redshift bin 
z— 0-0.2 for three mass bins of the parent clusters AAf : < 
9x 10"Mq, 9x 10" < M < 1.8 X 10^^ Mq, and 1.8 x 10^^ < 
Af < 3.6 X 10^^ Mq, which are consistent with the luminosity 
bins adopted to draw the observed statistics (Giovannini et 
al. 1999; Giovannini & Feretti 2002). 

First we run a large number, Af, of trees for differ- 
ent cluster masses at 2; = 0, ranging from ~ IO^'^A/q to 
~ lO^AfQ. Thus, for each M, we estimate the formation 
probability of radio halos in the mass bin AM as : 



rAM 
JM 



(44) 



where tu is the time that the cluster spends at z < 0.2 in 
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the mass bin AM with x ^ Xmin*^ and td is the time that 
the same cluster spends in AM with x < Xmin- 

Thus the total probability of halo formation in the 
mass bin AAf is obtained by combining all the contribu- 
tions (Eq. 44) weighted with the present day mass function 
of clusters. 

We consider two possible cluster mass functions: the 
Press & Schecther mass function (1974, P&S): 



nps{M,z)dM 



2 j}^ Scjz 



(M) 



da{M) 



dM 



exp 



2a2(M) 



dM . 



(45) 



were a^{M) is given by Eq.(5) and the other quantities are 
given in Sect. 3.1, and the Sheth & Tormen (1999, S&T) 
mass function which is obtained from a fit to numerical sim- 
ulations and which predicts smaller and larger values of the 
cluster number density for small and large masses, respec- 
tively. We checked that the probability to have a radio halo 
obtained by making use of the P&S and S&T mass func- 
tions are consistent within few percent for the considered 
mass bins. 

In Fig. 7 we plot the occurrence of radio halos with a 
typical radius Rh ~ 500 kpc, as a function of 774 compared 
with the observated statistics (see caption) . We find that the 
relatively high occurence of radio halos observed in massive 
clusters can be well reproduced by our modelling under very 
reasonable conditions, i.e. that a fraction of 20-30% of the 
energy of the turbulent motions (about few percent of the 
thermal energy) is in the form of compressible MS waves. 
In addition, we find that there is a range of values of the 
parameter r]t (0.2 < rft < 0.26, for B ~ 0.5/^G) for which the 
theoretical expectations are in agreement with the observed 
statistics in both the considered mass bins: ~ 30% and ~ 4% 
in the high and medium mass bins considered, respectively. 
Finally, we find that the expected probability to form giant 
radio halos in smaller clusters (not reported in Fig. 7) is 
negligible, in agreement with the observations. 



8 SUMMARY AND DISCUSSION 

Crucial constraints on the origin of radio halos are provided 
by statistical studies which show a connection between the 
formation of these sources and cluster mergers, and also find 
an abrupt increase of the occurence of radio halos with the 
mass of the parent clusters. 

- The first goal of the present paper is to check if cluster 
turbulence generated during mergers may be able to drive 
efficient particle acceleration processes in the ICM. 

- The second goal, in the framework of the turbulent- 
acceleration hypothesis, is to investigate if the hierarchical 
formation process of galaxy clusters can naturally account 
for the increase of the radio halos' occurrence with cluster 



"I Since clusters in our synthetic population never liave x >> 
Xmin, the condition x > Xmin guarantees a synchrotron spectral 
index compatible with that of radio halos. 
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Figure 7. Expected formation probability of radio halos {Rh — 
500 kpc, B ~ 0.5/iG) in a EdS cosmology as a function of pa- 
rameter rjt in two different mass bins (solid lines with error bars): 
binA= [1.8 - 3.6] 10^^ Mq and binB= [0.9 - 1.8] W^'^Mq. The 
two bottom dashed lines mark the observed probabilities for ra- 
dio halos in the mass binB while the two top dashed lines mark 
the observed probabilities in the mass binA. The two reported 
observational ranges account for Icr errors. The theoretical errors 
are estimated by extracting sub-samples of galaxy clusters from 
the synthetic population with a Monte Carlo procedure. 



To achieve these goals we developed a statistical 
magneto-turbulent model which is based on the following 
steps : 

i) Extensive merger trees of galaxy clusters with different 
present day masses are obtained. The trees are calculated 
making use of a procedure of Binary Merger Tree Method 
which is based on the extended PS formalism (Sect. 3). The 
temperature of the ICM is estimated at each redshift from 
the observed M-T relationships. 

li) Cluster turbulence is assumed to be injected during 
cluster mergers by the crossing of the infalling subclusters 
into the larger ones. To be conservative, turbulence is as- 
sumed to be injected in the major subcluster only within 
the volume sweeped by the minor subcluster (Sect. 4). The 
injection rate of MS waves is assumed to be a fraction, rjt, 
of the turbulence injection rate. 

The injection spectrum of MS waves is assumed to be 
a simple power law which extends over a broad range of 
scales (Sect. 4), or a delta-function from which turbulent 
cascade is originated (Appendix B). In both cases the max- 



imum/injection scale is fixed at Lt 



2rs. The resulting 



spectrum of MS waves is then calculated assuming station- 
ary conditions within a crossing time for each merger and by 
taking into account the relevant damping processes (or cas- 
cading processes) in the ICM (Sect. 4). The evolution of the 
spectrum of MS waves during cluster formation is calculated 
by combining the effect of all mergers. 

in) The evolution of relativistic electrons in galaxy clusters 
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is calculated considering the acceleration by MS waves (ti) 
and the energy losses. Relativistic electrons are assumed to 
be continuously injected in the ICM by shocks, AGNs and 
star forming galaxies in the clusters during their life. The 
total energy budget injected in the relativistic electrons is 
assumed to be a fraction, rje, of the thermal energy of the 
clusters at the present epoch. We do not follow the evolu- 
tion of the relativistic hadronic component since the most 
important damping of MS waves is with thermal electrons 
(Sect. 4.2) and thus the relativistic liadrons do not affect 
significantly the electron-acceleration process. 

To match the redshift range spanned by observational 
studies we calculate the model expectations for z < 0.2. The 
comparison between model and observations is performed in 
two main steps : 

i) First we consider the case of a typical massive clus- 
ter of our synthetic population and calculate the expected 
synchrotron and inverse Compton emission as a function of 
rjt and rje- We find that the typical radio luminosity of rar 
dio halos and the HXR luminosities can be obtained by our 
model provided that a fraction of the cluster turbulence, 
rjt ~ [0.2 - 0.3] (B^g/0.5)"^/^ (B^g being the volume aver- 
aged field strength within Rh in units of fiG) , is channeled 
into MS waves during major mergers and that the energy 
injected into relativistic electrons is 10~^ — 10~* times the 
present energy of the thermal pool (Sect. 6.2, see also the 
discussion below). 

n) Then, we compute the occurence of radio halos with 
the mass of the parent clusters. More specifically, we cal- 
culate the minimum particle acceleration coefficient, \mim 
which is required to efficiently boost the accelerated elec- 
tron population and produce radio emission with the spec- 
tral slope typical of radio halos. We thus identify the galaxy 
clusters containing a radio halo as those clusters in our syn- 
thetic population for which x > Xmin (see Sect. 7 for de- 
tails). The radio halos' occurrence is calculated in three 
mass bins consistent with those adopted in observational 
studies (< 9 • lO^Mgh^oS 9 • lO" - 1.8 • lO^^Moh^oS and 
1.8 ■ 10^''' - 3.6 ■ lO^^Mohso^). We find that for a single range 
of values of rit it is possible to account for the observed 
probabilities in all the three mass bins: about ~ 30% and 
~ 4% in the larger and medium mass bins, respectively, 
while the probability to find a radio halo in a cluster with 
mass < 9 • W^^MqIl^q is found to be negligible. 

As a general conclusion we find that the model expec- 
tations are in good agreement with the observational con- 
straints for reliable values of the two free model parameters; 
Vt, Ve- 

We also find that given these parameters and the phys- 
ical conditions in the ICM, the cascade time of the largest 
eddies of the MHD turbulence is of the order of ~1 Gyr. 
Consequently the difi'usiori and transport of these large scale 
eddies and waves may give a fairly uniform turbulent inten- 
sity within a relatively large volume (> Rh)- Finally, we find 
that the two extreme scenarios considered in our model, i.e. 
an injection of the MS waves with a power law spectrum, 
or with a single scale delta-function, provides very similar 
results since the process of particle acceleration basically de- 
pends on the energy flux injected into MS waves (which is 



dissipated at coUisionless scales) and on the physical condi- 
tions in the ICM (Appendix B). 

Thus, although the necessary approximations adopted 
in our formalism, we have shown that particle acceleration 
processes, which are invoked to explain the morphological 
and spectral properties of radio halos, can also account for 
the statistical properties of this class of objects. 
The following items need some further discussion: 

• An important finding of our calculations is that only 
massive clusters can host giant radio halos {Rh > 500 kpc) 
and that the probabilty to form these diffuse radio sources 
presents an abrupt increase for clusters with about M > 
2 X 10^^ Mrs,. 

Fig. 4 shows that the energy of the turbulence injected 
in galaxy clusters is expected to roughly scale with the ther- 
mal energy of the clusters. This seems a reasonable finding 
which immediately implies that the energy density of the 
turbulence is an increasing function of the mass of the clus- 
ters, £t ocT oc M". In addition, in the case of clusters with 
mass M < lO^^M© the infall of subclusters through the 
main one injects turbulence in a volume Vt smaller than that 
of giant radio halos, Vh, and thus the efficiency of the mech- 
anism is reduced by about a factor of Vt/VH (Sect. 4.1). On 
the other hand, major mergers between massive subclusters 
are expected to inject turbulence on larger volumes, of the 
order of Vh (or larger, e.g. Fig. 3b), and thus the efiiciency 
of the generation of radio halos is not reduced. 

More quantitatively, focussing for simplicity on what 
happens during a single merger event, the efficiency of the 
particle acceleration in the fixed volume Vh = 'inR^/S 
can be derived from Eq.(40) : x « 9{rs, Rh){M/ Rf^^ /VT, 
where the term T^^^'^ is due to the stronger damping of MS 
waves on thermal electrons with increasing the temperature 
of the ICM (Eq. 20). Thus the acceleration efiiciency within 
Vh is found to scale about with x M^~''^'^g{ra, Rh) oc 
^0.75-1.25 75 y 3.10^^ Me, 1.25 for M < 10^^' Mq). 

• Several mechanisms can provide injection of turbu- 
lence in the ICM during cluster mergers. We have just fol- 
lowed a simple approach which allows us to estimate the 
injection of turbulence during the crossing of smaller clus- 
ters through the more massive ones. It should be reminded 
that in the calculations we have adopted a typical radius of 
a radio halo, Rh ~ 500 kpc, and assumed that turbulence 
injected in a smaller volume is diffused on the scales of the 
radio halo, while the effect of the turbulence injected out- 
side Rh is not considered. However, the stripping radius, in 
the case of major mergers between very massive subclusters, 
can be larger than Rh ~ 500 kpc and thus the turbulence 
injected by these massive mergers can power particle accel- 
eration also on larger scales. The relativistic electrons ac- 
celerated at these scales can significantly contribute to the 
IC spectrum and thus the IC luminosities given in this pa- 
per (e.g.. Fig. 6) may be underestimated. On the other hand, 
the volume integrated synchrotron spectra should be mainly 
contributed by the emission produced within Rh due to the 
expected decrease of the magnetic field strength with radius. 
Thus our results, which are essentially based on the syn- 
chrotron properties of radio halos, should not be affected by 
the presence of non thermal emission from very large scales. 

• An important result of this work is that the en- 
ergy which is required to be injected in relativistic electrons 
in the volume of galaxy clusters is of the order of a few 
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10^'* X (_B^g/0.5)^^ of the present day thermal energy of the 
ICM. This value basically depends on the balance between 
the electrons' energy losses and the turbulent-acceleration 
efficiency which is experienced by the relativistic electrons 
injected in the ICM during the last few Gyrs. Since our 
calculations are performed by assuming the physical con- 
ditions of the ICM as averaged over the cluster volume, 
the required injected energy in relativistic electrons may 
be substantially higher in the central regions of the clus- 
ters where the high density of the thermal plasma causes 
stronger Coulomb losses. We notice that the required values 
of r]e can be easily provided by considering the injection of 
relativistic electrons in the ICM from AGNs, galactic winds, 
and large scale shocks (e.g., Bicrmann et al., 2003 for a re- 
view) . 

The requested values of the energy injected in relativis- 
tic electrons in the ICM are calculated thorugh the paper by 
assuming s — 2.5 and pmin/mc — 60 (Sect. 5.1). The results 
however should not be very sensitive to these assumptions, 
and they would be only sensitive to the total number of rela- 
tivistic electrons injected in the ICM during the cluster life. 
Indeed, the turbulence experienced in the ICM basically in- 
creases the cooling time of the injected electrons which are 
then mantained at the peak of their cooling-time curve (i.e., 
at 7 ~ 100 — 200, e.g. Sarazin 1999) and thus boosted at 
higher energies during an efficient re-acceleration period. In 
order to test the poor dependence of our results on the as- 
sumptions on s and pmin, wc re-calculate the value of rje 
by assuming s = 2.2 — 3.0 and pmin/mc = 20 — 100. We 
find that different assumptions requires values of rje within 
a factor of ~ 3 to reproduce a given synchrotron power. In 
particular we find that rje decreases with increasing s (or 
with decreasing pmin). 

It should be stressed that the amount of injection of rel- 
ativistic electrons required by our model is orders of mag- 
nitude smaller than that needed by models which assume 
a simple contirmous injection of a power law energy distri- 
bution of the emitting electrons in the ICM (e.g., Sarazin 
1999). This is mainly because during an efficient accelera- 
tion period the spectrum of the relativistic electrons is not 
a steep power law in which the bulk of the electrons is at 
low energies. During this period the bulk of the electrons, 
accumulated at 7 ~ 100 — 300, is boosted at higher energies 
and essentially piled up in the energy range responsible for 
the synchrotron emission in the radio band. 

• Since the present work is not aimed at reproducing in 
detail the properties of radio halos, in our calculations we as- 
sume that the magnetic field strength (within Rh) is roughly 
constant (B ~ 0.5/iG is assumed to constrain rjt and rje). 
Larger values of B (but still under the conditions in which 
the radiative cooling of electrons is dominated by IC emis- 
sion, i.e. B < 3/iG) would allow to radiate the synchrotron 
photons at higher fi-equencies (Sect. 6.3, Eq. 43) and this 
would imply that lower values of rjt (rjt oc B^^^^, Sect. 6.2) 
are required to form radio halos. On the other hand, the 
discovery of HXR tails in galaxy clusters has revealed that 
the non-thermal spectra of these objects are dominated by 
the IC component which has a luminosity ~ 10^ times larger 
than the synchrotron component. If confirmed, these obser- 
vations indicate that the volume-averaged magnetic field 
strength should be < 0.5/LtG (e.g., Fusco-Femiano et al., 
2003). However, as discussed above, a relevant contribution 



to the IC spectrum of galaxy clusters can be provided by 
electrons accelerated by turbulence injected in the outer re- 
gions (> Rh) and thus values B>l;uG in the synchrotron 
emitting volume may be still compatible with the observed 
IC components. 

• As already stated in the model calculations we have 
assumed a typical mean radius of the radio halos and a value 
of the magnetic field strength B which are independent from 
the mass of the parent clusters. If radio halos in more mas- 
sive clusters arc larger than those in smaller ones, then this 
approach should underproduce the expected probability to 
find radio halos in the smaller clusters with respect to the 
larger ones. The fact that the values of rjt required to match 
observations in the intermediate mass bin arc found to be 
slightly larger than those in the more massive bin (Figs. 7 
and Al) may reflect this effect. On the other hand, if B in- 
creases with the mass of the parent clusters (with B < S/iG) 
then the synchrotron emission would be boosted at higher 
frequencies and the expected probability to find radio halos 
in the case of larger clusters would be slightly increased with 
respect to our present expectations. 

• Finally, the model results obtained with a EdS cos- 
mology have been compared with those obtained assuming a 
ACDM cosmology (Appendix A) . Wo find that the possibil- 
ity to explain the observations in the redshift bin — 0.2 does 
not depend critically on the adopted cosmology. In particu- 
lar, although the model is found to be slightly less efficient 
in a ACDM cosmology, also in this case the occurence of 
radio halos can be matched for viable values of the param- 
eters, and a single range of r)t is found to be able to explain 
observations in all the mass bins; this would strengthen our 
conclusions. 

Future studies with radio (LOFAR, SKA) and hard X- 
ray (ASTR0-E2, NEXT) observatories will be crucial to 
constrain the occurrence and evolution of the observed non- 
thermal diffuse emission in galaxy clusters and thus to per- 
form a detailed comparison between observations and model 
expectations. 
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APPENDIX A: THE CASE OF ACDM MODEL 

In order to have a prompt comparison with observational 
studies, all the results given in the paper are obtained in 
a Einstein-De Sitter (EdS) cosmology. In this Appendix we 
adopt a ACDM cosmology (we use O,m{0) = 0.3, nA(0) = 
0.7, (78 = 0.9 and ho = 0.7) and discuss the differences with 
the EdS case. 

In the ACDM cosmology the critical overdensity (in 
Eq.l) as a function of the cosmic time is given by (Kitayama 
& Suto 1996): 



Sc{t) = { 1 + 0.0123lognm{z] 



(Al) 



where Qrn{z) is the mass density ratio at the redshift z. 



The growth factor in Eq.(Al) is (Peebles 1980, Eq.13.6): 



(A2) 



D{x) 



(^3 + 2)V2 r ^3/2(^3 ^2)-3/2^^ (A3) 

Jo 



„3/2 
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Figure Al. Expected formation probability of radio halos (Rh — 
as a function of parameter r]t in a EdS 
cosmology (solid lines with error bars) and in a ACDM cosmol- 
ogy (dotted lines) in the mass bins; binA=[1.8 — 3.6] 10^^ Mq K^^ 
and binB=[0.9 - 1.8] 10^^ Mq for EdS case and binA=[1.9 - 
3.8] • lO^S Mgh^Q^ and binB=[0.945 - 1.9] • lO^^ ^Q^^o for the 
ACDM model. The two bottom dashed lines mark the observed 
probabilities for radio halos in the mass binB while the two top 
dashed lines mark the observed probabilities in the mass binA; 
observational regions account for la errors. 



where xq = {^Y^^ and x = xo/{l + z). 
The ratio of the average density of the cluster to the mean 
density of the universe at a given z, Ac{z) (in Eq.(lO)), in 
the ACDM model is given by (Kitayama & Suto 1996): 



(A4) 



where uj{z) = ^lm{z)~^ — 1. 

Following the procedures adopted in the case of the EdS 
cosmology, we compute merger trees (Sect. 3), turbulence in- 
jection rate and spectrum of the MS waves (Sect. 4), particle 
evolution (Sect. 5) and non-thermal emission (Sect. 6) from 
galaxy clusters and thus the expected formation probability 
of radio halos for z < 0.2 (Sect. 7). In Fig.Al, we report 
the comparison between the probability to form radio halos 
obtained in the two cosmologies. The comparison is derived 
by converting the virial mass of the clusters from a EdS into 
a ACDM model: 



M 



[Ac{t)pm{t)]EdS 
[Ac(f)pm(t)]A 



1/2 



(A5) 



where pm is the mean mass density (dark and barionic) of 
the Universe. Thus the calculations with a ACDM model are 
performed for the mass bins [0.945 - 1.9] • 10^^^ Mohfo^ and 
[1.9-3.8] • 10^^ M0hfo\ 

As expected, we find that at 2 < 0.2 the results are rela- 
tively independent from the considered cosmology, with the 
ACDM model being only slightly less efficient. In particular, 
as in the EdS case we note that it is possible to find a unique 



interval in r/t in which the model reproduces the observed 
halo formation probability for both the cluster-mass bins. 

In the ACDM Universe the structures start to grow at 
early time with respect to the EdS case, the merging rate at 
z < 0.2 is consequently reduced, and thus particle acceler- 
ation is less efficient. However, this is roughly compensated 
by the fact that in a ACDM Universe the observed radio 
halos are smaller and less luminous than in the EdS case. 



APPENDIX B: TURBULENCE INJECTION AT 
A SINGLE SCALE 

In this Appendix we adopt the scenario in which MHD 
turbulence is injected in the ICM at a large single-scale, 
kmin ~ 7r/'"s, fiom which the MHD turbulence cascade is 
originated. 

The mean free path, Lmfp, in the ICM marks the 
boundary between the coUisionless regime {k > 2Tr/ Lmfp) 
and the coUisional regime (fc < 2n/Lmfp)- It is given by 
(e.g., Braginskii 1965): 



(kpc) - 



300(^)^(^- 



(Bl) 



which, for the typical values of the cluster temperatures and 
of the mean thermal density within Vh, is of the order of 
100 kpc. 

Once MS waves are injected at kmin, the process of 
wave-wave coupling generates a turbulence cascade. The 
cascade time of fast MS waves at the wavenumber k is given 
by (e.g., Yan & Lazarian 2004): 



rkk{k) 



vm 
kvl 



VMp 



(B2) 



k^k 

so that the diffusion coefficient in Eq.(19) is given by : 

k^ v\ k^W, 

Tfcfe VM 2Wb 

In the quasi linear regime, the spectrum of the waves due to 
the cascading process can be calculated solving Eq.(19) and 
neglecting the contribution due to the damping terms : 



dt 



dk 



(B4) 



with /k = IoS{k — kmin) and h ^ rjtvf p^Tvrl /Vh) (Sect. 4.1). 
The stady state solution of Eq. (B4) is a Kraichnan-like spec- 
trum : 

'2p/„«AfNl/2 



^2pioVM_y 



-3/2 



(B5) 



This spectrum extends down to a truncation scale at which 
the cascading time, r^k, becomes substantially larger (i.e., ^ 
times, ^ ~ 1 — 3) than the damping time scale, Td ~ ^t^he 
(Eq. 20). In the coUisionless regime, this truncation scale, 
Ltr ~ 2ii/ktr, is obtained from Eqs.(20), (B3), and (B5), 
one has : 



0.23 r , 



lO^km/s' 



4_R| 
3 ri 



(B6) 



which typically falls in the range 10-30 kpc for our synthetic 
clusters (note that such scale is smaller than Lmfp and thus 
the estimate can be done under the assumption of a col- 
lisionless regime), i.e. a factor of 30-100 smaller than the 
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value of the typical turbulence injection scale (Sect. 4.1-4.2, 
Fig. 3b). 

The picture of the model in this Appendix is thus that 
the injection of MHD turbulence occurs at a maximum scale 
of the order of 1 Mpc which is larger but relatively close 
to the scales typical of the coUisionless regime. The wave- 
wave coupling then leads to a power-law inertial range with 
a Kraichnan spectrum which is approximatively mantained 
down to ~ 10—30 kpc. At these scales the damping time with 
the thermal electrons becomes considerably shorter than the 
cascading time-scale and the turbulence cascade is broken. 

Under these conditions the acceleration time of rela- 
tivistic electrons, race, is dominated by the contribution from 
the spectrum of the waves at the truncation scale and it can 
be obtained from Eqs.(35) and (36) : 



X oc VMWk{k = ktr)kl 



(B7) 



An important point is to check if the scenario adopted 
in Sect. 4 and that adopted in this Appendix give consistent 
results. In the scenario adopted in Sect. 4.2, the spectrum 
of the MS waves is approximately given by : 



(B8) 



and thus, since lo = J hdk and the damping time scale is 
Td{k) ~ r~j^ oc k~^, one has : 



Wkkdk ~ IoTd{ktr)ktr 



(B9) 



APPENDIX C: SYNCHROTRON AND IC 
RADIATION 

In this Appendix we briefly describe the formalism adopted 
to calculate the synchrotron and IC emission from galaxy 
clusters. 

The synchrotron emissivity is given by: 



"^^y [ f ' dpd9sin''9N{p,t)F{—) (CI) 
Jp Jo 

where F(^) is the synchrotron Kernel given by (e.g., Ry- 
bicki & Lightman, 1979): 



K5 {y)dy 



(C2) 



with i^c = {3/4:'K)p^eBsinO/{mc)^ 



The IC emissivity due to the scattering off the CMB 
photons by relativistic electrons is given by (e.g., Blumen- 
thal & Gould, 1970): 



Jicivijt) = {2nrom'^c)^hi'i / / dvdp 



N{p,t)p- 



where = 2.73(1 + x) is the temperature of the CMB 
photons. 



where Cpp does not depend on the turbulence spectrum and 
energy (Eq. 35). On the other hand, in the scenario adopted 
in this Appendix the spectrum of the MS waves is approxi- 
mately given by : 



(BIO) 



and thus, since the cascading time scale is Tkk oc k one 
has : 



Wkkdk ~ 2/o T/jfc(/u — ktr^^tr 



(Bll) 



as a consequence, since ktr is the scale at which the damp- 
ing time scale and the cascading time scale are comparable, 
we expect that the two scenarios would provide a similar 
acceleration efficiency. 

More specifically, we can calculate the electron acceleration 
coefficient due to a single merger event (with > Rh) in 
the framework of the scenario adopted in this Appendix. 
Prom Eqs.(B5), (B6), (B7), (9), and the expression for 
given in this Appendix, one finds : 



X~2xl0-'''^??t 



M 



2 X 1015 M( 



3/2 



\7keVj 



-1/2 



{rs/500kpcf 



(ii|/500fcpc) 



(B12) 



which is close to the value given in Eq.(40), and thus proves 
that the main results of our model do not crucially depend 
on the assumptions on the specific injection process (and 
spectrum) of the MS waves. 



© 0000 RAS, MNRAS 000, 000-000 



